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Abstract 

A systematic study of a wide class of nonrelativistic models of 
bb quarkonia is described. It is found that the potential V(r) = 
0.706380(V^- a46 r 0442 ) + 8.8 1 71 5 (all in GeV) with the 5-quark mass 
nib = 4.80303 GeV gives a satisfactory description of the experimental 
data below the threshold for strong decays (x 2 /DF = 6.5/7). Limi- 
tations and implications of this observation are discussed. 
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Dozens of nonrelativistic and relativistically improved models of bb quarko- 
nia have been published. Many references can be found in the recent review 
[|]. A fully relativistic treatment is beyond our reach, even for the much sim- 
pler e + e~ system, but the folklore is that good results can be obtained with 
a nonrelativistic potential, if the hamiltonian is interpreted as an effective 
hamiltonian with the coefficients renormalized by the relativistic corrections. 
Explicit relativistic corrections are necessary to describe some purely rela- 
tivistic effects, like the fine splitting. If, however, one limits the discussion 
to centres of gravity of the multiplets, than at the phenomenological level, 
as seen e.g. from [|TJ, the effective nonrelativistic models are doing about as 
well as the relativistically improved ones. In the literature a model is con- 
sidered good, if it reproduces e.g. the masses of the quarkonia within a few 
MeV. The experimental errors on the masses today [Q], however, are of the 
order of 0.2 MeV. This rises the interesting question: is it possible to fit the 
data quantitatively, i.e. within the experimental errors and if so, what is the 
corresponding effective potential. The answer seems useful for at least the 
following two reasons. Quark - antiquark potentials derived from studies of 
quarkonia are used in a variety of applications. Let us mention as examples 
applications to heavy light systems , to be mesons || , and to it production 
M. It is obviously advisable to use as good potentials as possible and to 
know what are their limitations. On the other hand, if it is not possible to 
fit the data with a nonrelativistic model, this may be an interesting hint on 
how to construct a relativistic theory. 

The first difficulty is that there is no such thing as a standard nonrela- 
tivistic quarkonium model. The Schrodinger equation is, of course, 

- — V 2 ^(f) + V(rU(r) = Ei/;(r), (1) 

but the quark mass m;, and the potential V(r) vary from paper to paper. 
Our strategy is to limit the discussion to potentials of the form 

V(r) = -ar~ a + br? + Ct, (2) 

where a, b, a, (3, Ct are nonnegative constants. At least ten potentials of this 
general form, but with various values of the parameters, have been proposed 
in the literature. Thus the famous Cornell potential || has a — (3 — 1. The 
potential advocated by Lichtenberg and collaborators [|7| has a = (3 = 0.75. 
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The potential used by Song and Lin in ref. || has a = (3 = 0.5, while Song's 
potential used in ref |J has a = (3 = |. The logarithmic potential of Quigg 
and Rosner [|1(J corresponds to a = (3 — > 0. Potentials with a ^ (3 have also 



been popular. Thus Martin JTTJ] has suggested a = 0, (3 = 0.1, while Grant, 



Rosner and Rynes |T2| prefer a = 0.045, (3 = 0. Heikkila, Tornquist and 



Ono [13] tried a = 1, (3 = Some very successful potentials known from the 



literature are not of this type. Examples are the Indiana potential [TJJ and 



the Richardson potential |15| . Our first observation is that these potentials 
are very similar to the potentials considered by us, if their free parameters 
are fitted to the data as explained below. This is illustrated in Fig. 1, where 
the Richardson potential and the Indiana potential are compared with two 
potentials of type (0). 

Let us make three comments concerning this figure. A figure analogous 
to our Fig. 1 has been published long ago by Buchmiiller and Tye jTJ 



The agreement between the curves on their plot was not as good as in Fig. 
1, in spite of the fact that they adjusted the constants in the potentials 
so as to make all the potentials coincide in a chosen reference point. The 
reason is that we have been much more selective, than was possible at their 
time, in the choice of the "good" potentials for comparison. Our result 
supports the conjecture of Quigg and Rosner |DJ, who concluded from their 
inverse scattering analysis that the first few L = energy levels and the 
corresponding leptonic widths determine to a large extent the potential in 
the region relevant for the quarkonium calculations. In our fits we used 
more observables, but very similar results can be obtained by using just the 
masses. Thus, a version of our analysis is similar to that from ref. [|17] 



except that instead of the leptonic widths we are using the centres of gravity 
of the L = 1 states. Note, however, that we have no proof that a completely 
different potential would not fit the data as well. Finally, let us stress that in 
the regions of very large and very small values of r, not show in the figure, the 
potentials are completely different from each other, but this has practically 
no effect on the calculations for the quarkonia. 

It is convenient to rewrite equation (0) in the reduced form 

-V 2 0(p1 + V(p)0(p)=^0(p), (3) 

where 
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P = ^, (4) 

^ = (a) ' ® 

V(p) = C(pP-p- a ), (6) 

£ = ^(E-Ct), (7) 

C = m b a\ a - 2 . (8) 

The idea is to concentrate on the observables, which depend on the parameter 
C only. We choose 



t _ M(3S)-M(2P)_ a774±a006 _ 



M(25) - M(1P) 
Af(25) - M(LP) 
M(2S) - M(1S) 
lfe(0)| 2 



0.219 ±0.001, (11) 



\M0)\ 2 

1^35(0) | 2 



0.492 ±0.111, (12) 
0.433 ±0.071, (13) 



l<M0)| 2 

fee = |Vis(0)|s<lP|r|25) = 2.29±0.16, (14) 
b 7 = \^ 1S (0)\I (2P\r\3S) = 1.59 ±0.15, (15) 

" 8 = |^iy = ' 110±0 ' 009 ' (16) 

Thus the x 2 distribution corresponds to seven degrees of freedom. All the 
numerical values are calculated from the data given in the 1994 Particle Data 
Group Tables 0. 

This choice of observables requires some comments. Since our model is 
nonrelativistic, we have replaced the masses of the \ states by the centres of 
gravity of the multiplets. It would have been nice to be able to include in the 
averagings also the masses of the rjb and the 1 Pi states, but these masses are 
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unknown. We have ignored the quarkonia with masses above the threshold 
for strong decays. Their analysis would require coupled channel calculations, 
which are much less well-defined. The extraction of values of the wave func- 
tions at the origin |^(0)| from the experimentally measured leptonic decay 
widths is based on the Van Royen - Weisskopf formula with a first order 
radiative correction (cf. e.g. Since the first order correction is about 

30%, the second order correction is probably significant. Unfortunately it is 
not known. Moreover, the formula itself leads in certain cases to paradoxes 
(cf. jn| and references given there). It is believed that the resulting uncer- 
tainties largely cancel in the ratios 64 and 65. In the observables be and 67, 
however, they introduce a systematic error of perhaps some 7%. This has not 
been included in our quoted errors. Thus at this point we underestimate the 
confidence levels. Finally, the relation between the dipole matrix elements 
occurring in the observables b 6 ,b 7 ,b 8 and the measured dipole transitions is 
for quarkonia less close than for atoms (cf. e.g. [[0|], pop ). For our best fit, 



and for the best fits to some other potentials, we have recalculated the full 
matrix elements without the multipole expansion. This yields the correction 
factors, which for transitions between S and P states are 



r , l <PM k i)i\s>? + 2\<vM k i)i\s>\- 

p " l<pl>>l ' ( } 

where s denotes the S"-wave function, p denotes the radial part of the P-wave 
function, ji are the spherical Bessel functions and k is the length of the wave 
vector of the emitted photon. The corrections for the transitions 2S — > IP 
and 3S — > 2P related to b 6 and b 7 are below 0.5%, i.e. negligible compared 
to the experimental uncertainties. The ratio of the transition probabilities 
(2P — > 1S)/{2P — > 25*) related to b 8 increases significantly, making our fits 
worse, but not bad. The x 2 of the overall fit increases by about two units. 

As a first step we have calculated the observables &i , . . . , b 8 and the corre- 
sponding values of x 2 f° r many of the existing models. Since our comparison 
with experiment eliminates three of the four parameters and adjusts the 
fourth to fit the data as well as possible, our agreement with experiment is 
usually better than in the original papers. A representative selection of the 
results is given in Table 1. We conclude that none of the models known to 
us from the literature fits the data in the sense of the \ 2 test. Therefore, we 
have explored the quality of the fit in the region < ot < 1.2, < /3 < 1.1 
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of the a, (3 plane. The resulting map is shown in Fig 2. We have plotted the 
parameter log Thus the large variability of \ 2 is clearly visible. We have 
also marked the points corresponding to some of the models known from 
the literature. Good fits correspond to the region \ 2 < 7 i.e. to the region 
shaded on the plot. All the models miss this region, though some come fairly 
close to it. 

In order to get a good fit we choose in the shaded region a point (see 
Figure 2) with the simple coordinates a = 1, (3 = 0.5. Minimizing \ 2 with 
respect to the constant C and then choosing the constants A, mj, and Ct so 
as to reproduce correctly the leptonic width of the T(1S) state, i.e. |?/>is(0)| 2 , 
the mass difference M(3S) — M (15), and the mass M(1S) and we find after 
substitutions 

V{r) = 0.706380 L/F - j + 8.81715, (18) 

where V(r) and r -1 are in GeV. The constants are given with the precision of 
six digits in order to assist the reader, who would like to check our calculation. 
The corresponding quark mass 

m b = 4.80303 GeV (19) 

is quite reasonable. The corresponding predictions for the parameters b±, . . . , bi 
are shown in the third line of the table. As expected the fit is very good - 
corresponding to a confidence level of 48%. A more precise treatment of the 
parameter b 8 , as mentioned above, yield x 2 = 8-8 which corresponds to a 
confidence level of 27%. The next best potential, the Indiana potential []Tj], 
gets a similar correction, thus it does not become competitive. 

For r — > our potential has the r _1 dependence corresponding to one 
gluon exchange. With present data, however, we have no evidence for the 
additional factor l/log(Ar), which according to QCD should be introduced 
by the running of the coupling constant. The expected part of the potential 
linear in r is not seen. Probably the bottomonia are too small to reach 
sufficiently far into the asymptotic region of linear confinement. Perhaps a 
more flexible potential would exhibit the linear part. 

Our conclusion is that it is possible to reproduce with a nonrelativistic 
theory the observables b\ . . . , 6g within the experimental errors. The corre- 
sponding potential ([Tj|) and the corresponding estimate of the mass of the 
6-quark ( |T9"D are very reasonable. 
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Table 1 

Comparison of the predictions of some models with 

experiment. 





h 


b 2 


h 


&4 


h 


K 


b 7 


bs 


X 2 /DF 


PDG data 


0.6290 


0.7738 


0.2187 


0.49 


0.43 


2.31 


1.59 


0.110 




error 


0.0005 


0.0057 


0.0009 


0.11 


0.07 


0.16 


0.15 


0.009 




This paper 


0.6293 


0.7744 


0.2191 


0.49 


0.36 


2.26 


1.37 


0.124 


6.5/7 


Indiana 


0.6299 


0.7829 


0.2202 


0.48 


0.36 


2.18 


1.33 


0.124 


15.5/7 


Lichtenberg 


0.6283 


0.7950 


0.2172 


0.48 


0.36 


2.19 


1.34 


0.126 


26.3/7 


Richardson 


0.6276 


0.8106 


0.2150 


0.47 


0.36 


2.18 


1.34 


0.127 


74/7 


Song-Lin 


0.6382 


0.7246 


0.2375 


0.49 


0.35 


2.04 


1.21 


0.112 


850/7 


Cornell 


0.6128 


0.8951 


0.1946 


0.47 


0.37 


2.42 


1.55 


0.142 


2220/7 


Martin 


0.6363 


0.6891 


0.2707 


0.54 


0.38 


1.81 


1.06 


0.032 


3720/7 



Figure captions 

Figure 1 Dependence of four typical potentials on the distance r. The 
continuous curve correspond to potential (|18D ; the dotted curve corresponds 
to the Indiana potential flj§ and to the potential @ with a — (3 = 0.75 [[?]], 
which coincide at the scale of the figure; the dashed curve corresponds to the 
Richardson potential fll5 |. All the potentials have been scaled and shifted 
(see text) to fit the data. 

2 

Figure 2 Map of the parameter log y- in the a,P plane. The region of 
very good fits (x 2 < 7) is shaded. The black points correspond to: a = (3 = 1 
|], a = /3 = 0.75 i, a = (3 = § !, a = 1, (3 = § @, a = (3 = 0.5 §, 
TTfl, a = 0.045, = [HI, a = (3 -> (TO] and to 



= 0, (3 
1, (3 = 



0.1 



= 0.045, ^ = 0@, 
0.5 as proposed in the present paper. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503420vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503420vl 



